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Abstract 


The performance of a turbocompressor comprised of a large number of mech- 
anically independent compressor and turbine stages is considered. Examination 
of the equilibrium conditions of any stage shows its performance, for a given 
design, to be dependent on four variables and hence the performance of the com- 
plete unit is also dependent on the same four variables. Exact arithmetical 
calculation would be possible, but very tedious. In order to render computation 
easier, the conception of reflected conditions is introduced. Physically this implies 
independence of each stage, in some respects, from its neighbors. The assumption 
is likely to be in error least for intermediate stages, and most for the first and last 
stages. Methods of expressing the aerodynamic performance of the blades of the 
compressor and turbine stages are introduced to render calculation easier. Some 
calculated performance figures are compared with experiméntal values and show 
by their reasonable agreement that the assumptions used are valid. The methods 
can be applied to other types of complex turbocompressor units by suitable modi- 
fication of procedure and details. 


Introduction 


In 1929, A. A. Griffith put forward proposals for the use of the internal 
combustion turbine for aircraft propulsion (1). He showed that with the then 
existing knowledge such an engine was a practical proposition, and he also 
investigated mathematically the best arrangement of engine. A gas turbine 
engine consists essentially of an air compressor, combustion chamber, turbine 
to drive the compressor, and turbine to supply the useful output. Note that 
the two turbines may in fact form a single unit, and that in the common jet 
propulsion engine the second turbine is replaced by a propulsive jet. Griffith 
proved that in order to obtain a wide range of efficient operation it was desirable 
to arrange the compressor, and the turbine driving it, in a series of mechanically 
independent stages, each small compressor stage being driven by its own small 
turbine stage. The use of ‘small’ here denotes a stage over which the pressure 
ratio is small. The mechanical arrangement of such a unit isshown diagram- 
atically in Fig. 1. 

In this paper we discuss some aspects of the performance of such a unit, 
known as a contraflow compressor. The problem may be put as follows. 
Suppose we have such a unit of given geometrical design. We can express its 
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over-all performance in various terms, for example, the mass flow, and initial 
density of the air entering the compressor, and its compression ratio and effi- 
ciency, and the four corresponding factors for the turbine. We shall show that 
in fact there are but four independent variables, and they may be conveniently 
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Fic. 1. Diagram of a contraflow turbocompressor. 


taken as the mass flow and initial density of the air at the compressor inlet 
and at the turbine inlet. The problem then is to calculate from these four 
independent variables the values of the other four variables and also, the 
working conditions of any of the, possibly many, compressor or turbine stages. 


Nomenclature 


mass flow, pounds per second, 
annulus area, square feet, 
total pressure, poundals per square foot, 
static pressure, poundals ‘ 
total temperature, °K., 
static temperature, °K., 
density, pounds per cubic foot, 
peripheral velocity of blade, feet per second, 
gas velocity, 
gas axial velocity, 
gas whirl velocity, 
u/Va, 
mean blade radius, 
angle between vector representing a gas velocity relative 
to a blade and the axial direction, 
angle between vector representing a gas velocity relative 
to the frame of the machine and the axial direction, 
Suffix 1 indicates condition at entry to a stage, 
2 indicates condition at exit from a stage, 
Prefix c indicates a compressor stage, 
t indicates a turbine stage, 
n denotes conditions at stage number 7. 
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The stages are numbered from the H.P. end, i.e., the entry to the turbine 
and the exit from the compressor. 


th 


E.g., ntX2 indicates the value of quality X at exit from the ™ turbine stage. 


A bar over a quality denotes the value of the ratio. 
i = _ cd 
E.e.; xX aearest 

General Theory 


We start by considering the most general case. The assumptions we make 
are, that at any stage, we may typify air conditions over the annulus by refer- 
ring to conditions at the mean radius, and that in any stage there is negligible 
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Fic. 2. Contraflow compressor—typical velocity vector triangles. 


change in the axial velocity. This does not preclude a change over the whole 
compressor or turbine, but merely that such changes may be represented by 
small discontinuities between stages. 


In Fig. 2 we show diagrams representing the velocity vectors for any com- 
pressor stage and for any turbine stage. Then in any compressor stage n, the 
work done per pound per second of air flowing is, by Euler’s momentum theory 

nclencAY,, = ncllencVa [tanncar — tanncag], 


where Av,, is the change of whirl velocity. 


nc 
But tann.a1 <= + tann-B1 _ nceP + tann-B1, 


é 
nceVa 


and tanncBi= tancn+1)cB2=(n41pe@ — taNcn41ycQe. 








288 CANADIAN JOURNAL OF RESEARCH. VOL. 87, SEC. F. 


So the work done per pound per second = 
ncV"a-nc® [ne® + (n+1)cP — taNncae — tan¢n+1)ca2] . (1) 


Similarly in any turbine stage m, the work done per pound of gas per second 


is ntVa-- ntO[tannrae— tannray] . 
But tanniar = nid — tannrB1 
and tann1 = tangn—1 82 = tanin—1y1a2 — (n—1 td. 


So the work done per pound-second = 
ntV"a-nt@[tann—1) 1a2+ tannta2— (n—1 th — nt] . (2) 


For the rotor to be in equilibrium, the total work done by the turbine stage 
must equal that absorbed by the compressor, plus any frictional or windage 
loss. We may write the equilibrium condition as 

cM nd anceblnced + (nt1eb — tannca2— tann41) cao] = 
k.iM ni a-ntoltan(n—1ta2+ tanniae— (n—1t — nt). (3) 


Where .M and ;M are the compressor and turbine mass flows and k is a factor 
nearly equal to 1, representing friction and windage losses. 


Rearranging Equation (3) we have 


M. ae". n p [ 


k ne® + (n+1)cP — tannca2— tan(n+1) ca] > 


[tan(n—1)1a2 + taNnrae — (n—1) to — nt] 


i Sete Jaen ] M.¥.?. n 
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[tanncas + tanyn+1) cae]. 


This equation is the general equation relating the speed of any rotor to the 
speeds of rotors on either side, and can be used for any rotor. It is a second 
order difference equation and in its general solution has two arbitrary con- 
stants. It is evident that these two constants are related to the direction of 
the air before entering the compressor and the turbine, i.e., the effect of any 
stationary guide vanes. To solve Equation (4) as it stands, is very tedious, 
and for the purposes of this analysis, we can make certain simplifying assump- 
tions, the validity of which will be discussed later. 


Let us suppose 


we * De for all n, 
nd = for all n, 
nceao= C forl<2a<N, 
ntao= ft forl<n<wN 


Note that since @ = u/Vq = 7/0, then 0,°6 = Fda. 
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Equation (4) then becomes 
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then (n—1)cP . ke + neG|R1 + ko] + (n+1)c - k= ky . 2c + 2t. 


Ri. 2¢ + 2t 
Qkit ke] - 


Put ne@ = Ax" + 


Then substituting in Equation (5) we have 
Ax™""[Ro+(Ritks)x ot. kyx?] = (). 


So either Ax"—!= 0 which means that »-¢ is constant for all m, irrespective 
of any boundary conditions, which is obviously absurd, or else 


Rix? + (Rit ke)x + ko= 0 
x?-+- (1 + Re ky)x + ko/ki= 0. 
rr > - a ke 
rhe roots of this equation are x;}= —1, x2»= — = 

1 
So the general solution of Equation (5) becomes 
ky . 2c + 2t 


neo = A[—ko/ki]"+ B[—1]"+ 
¢ [ i] [—1] Q(k.t be) 


(6) 


It remains to determine A and B. We do this by considering rotor m = 1, 
and rotor” = N. Atn = 1, i.e., the first turbine stage, the direction of the 
air, relative to the frame of the machine, is given by tan™';, say. 

Then Equation (5) becomes 

eG Rit Ro] t+oc6. Ri= ki. 2c +t + Br. (7) 

At n = N, if the direction of the air, entering the compressor, relative to 

the frame of the machine is given by tan~'6,, then we have 


(n—1) cP . Rot neb(Rit Ro) = (c - B.)Rit 2t. (8) 
Substitution of Equation (6) in Equations (7) and (8) yields 
: .2c + 2 2\" .2c + 2 
(i+ k)| - he y aD Ri. ete 4 «| (®) A+B biter) 
ky 2(ki + ko) ky 2(Rki + Ro) 
= k\2c +t + Bt 


N-1 be 
ks |4 (- Bs) + B(—1)" tee > 
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h\N o . bic® +3) 
b+ BAP on k¥ oT | = (c — Bk + 2t. 
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Solving for A and B, and noting that, for an aircraft engine, N would desirably 
be even, so as to balance gyroscopic loadings, we have 
2kic + kiB-+ t os ki2c +t + B: 
ky ko 
1 + (k2/ki)* (2 + ko/k:) 
B _! — ke — Bil + R2/k:) ue 
ko(1 + k2/k:) 


These then are the two arbitrary constants to go in Equation (6). We note 
that Equation (6) gives us nc¢ as a function of 








n, the stage position and N the total number of stages, 
6; t, the gas leaving angles from the moving blades, 
Be, Bt, the gas entry angle to the compressor and turbine, 
ki, ko, i.e., sa vet and si 

r 


If we consider a given machine, c, t, B., 8:, 7, and k are all fixed. If now we 
specify M, and 0,, or what is the same, M and since 12 = M/pA Equation (6) 
is completely defined, and the value of .¢, and thence ¢¢ can be calculated for 
each stage of the compressor and turbine. Now for any given value of .M and 
the density at entry to the compressor, we can calculate y,v,, and since we 
know y-¢, y-¢@ is known. Once the blade velocity and axial flow velocity are 
known, the performance of the stage is completely fixed. Thus the conditions 
leaving the stage can be calculated and the next stage considered, and so on. 
In short, if we know the four independent variables M, p, and .M, .p the per- 
formance of the whole unit is completely specified. There is no loss of gener- 
ality in the simplifying assumptions made and this conclusion applies to any 
machine of this type. 


Reflected Conditions 


It is evident that detailed calculations are tedious, and, in any event, the 
blade angles in any real compressor will not be equal for all compressor stages, 
so that the mathematical expression for the parameter .,@ becomes most com- 
plicated. Fortunately a simpler method can be used for most cases, without 
introducing errors greater than those implicit in necessarily assumed data 
concerning blade performance, to which we will refer later. 


When n-81= ncB2, and niBi= niB2, we have a symmetrical or reflected vari- 
ation of the absolute velocity of the gas in its passage through the machine, 
which can be regarded as the condition to be aimed at in design. If we assume 
that this condition is satisfied, the mathematics simplifies considerably. Phy- 
sically, we are in effect considering each rotor on its own without influence 
being exerted on it by its neighbors. The justification for this rather drastic 
assumption is that it makes the sum practicable, and experience shows that 
the results given are sufficiently near the truth to be of great use. 
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Referring to the vector triangles of Fig. 2, it is evident that if we suppose 
8:= B2 in both cases, then the change of tangential velocity component can 
be written 


Av, = 2(.u — Uv, tan eae) 
tAv,, = 2( Va tan ;a2— 1) 


Expressing the equilibrium of the rotor, we have 








cM. .u.2(.u —Vq tan -a2) = k..Mu.2(w tan pa2— uM) 
cM.2.cbVa?(ed — tan ae) = k..M.2:d,,?(tan sa2— td) 
Vf a w2 
M$ v0 (.@ — tan Ao) = tan (a2 — oo 
k 
ae aS 
ages . tan cao+ tan ;a2 
or 2? = : 
Mov. , 1 
k ¢ 
Remembering that @ = t/ta = t/a 
and Va = M/pA, 
we have 
— . tan ae + tan ;@2 
ae ae m7 
M*r ] 
Me is 
kRAp rAp 
Hence for any given stage .¢ is a function of M and p, just as we found in the 
more general case, Equation (6). Also since ;¢ = = cP 
p eat 
tan ;a2 + ne tan .ae 
to = —____—-— ‘ (10) 
(+) 
k 


Note further that from the vector triangles and the assumed symmetrical 
conditions 


tan ,.a,;= 2. — tan ,a@o, (11) 
tan ,a,;= 2: — tan ;a2. (12) 
Finally we define wW = 2,6(tan sa2— 1d) 
and = 2.6(-¢ — tan <a). (13) 
The stage total temperature rise is AT = YBa (14) 
cKkp 


Ww . Va" 


tp 


and the turbine temperature drop is AT = (15) 
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Blade Performance 


It remains now to consider, in more detail, the aerodynamic performance of 
the blades. For this we rely on cascade tests (2) in which a pack of blades is 
mounted in a wind tunnel and the angle az, and the losses of pressure are 
measured for a range of values of a; and M,,;. Such tests are two dimensional 
and when the results are applied to predict the flow in an annulus some care 
is necessary. It is found by experience (3), that the values of the leaving gas 
angle measured in two dimensional tests may safely be applied to the annular 
flow, whilst to the two dimensional loss measured must be added due allowance 
for annulus and wall friction loss, and also an allowance for the fact that the 




















Fic. 3. Compressor blade performance curves—section Gottingen 436, 50 deg. Camber,—20.5 
deg. stagger. 


flow is in fact three dimensional and that there are slight variations in incidence 
up and down the blade, leading to secondary losses. These may be allowed for 
and in Fig. 3 we show a typical set of curves, relating the tangent of az to tan ai, 
and the Mach number. Fig. 4 shows the stage efficiency as a function of a; 
and Mach number. 


A third factor to be considered is the value of the axial velocity. This may 
be calculated from the mass flow and area and density to give a mean value, 
but in practice the growth of the boundary layer in an adverse pressure gradient 
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increases the axial velocity above this value. This may be allowed for if 
necessary by putting in an empirical factor multiplying the right hand side of 
Equations (14) and (15), known as a work factor, or it may be taken care of 
by putting in a factor multiplying the nominal flow areas by a so-called con- 
traction factor. Either, since the factors must be determined from analysis 
of test results, will yield reasonable results afterwards. The contraction factor 
is probably sounder physically. 
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Fic. 4. Estimated compressor stage efficiencies, assuming annulus and secondary losses to be 
150% of the minimum profile drag. 


Performance Integration 


We are now in a position to combine Equations (9) to (15) with the data of 
Figs. 3 and 4 and the corresponding data for the turbine blades to produce a 


complete performance analysis. 


The combination of Equation (11), with Fig. 3, yields a relation between 
e@, tan .a2 and ,.Mn, which is shown on Fig. 5. There is a corresponding rela- 
tion between ;@, tan ;a2, and ;Mn,, but in this case we do not show on Fig. 6 
any variation due to ;Mm, as in turbines of this kind, the turbine Mach 
numbers are always a long way from the critical values. From Figs. 5 and 6 
we can solve Equation (9). For typical values of A and 7, and M = 1.0 this 
is plotted on Fig. 7 as a curve of .@ against p. 


This is the key to the solution and there is no need to proceed in any detail 
further. 


For any contraflow compressor of specified geometry we can prepare curves 
such as Figs. 7 and 8 for each rotor. Then our four independent variables are 
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cM, .p, p, and M. The @ is most conveniently expressed at the low pressure 
end of the unit. Then since pj fixes .¢ and .M and ¢p fix ,vq, all the conditions 
at the L.P. rotor are known, in particular all gas velocities and directions. 


The value of .W gives us the temperature rise in the compressor stage, and 
Fig. 4 yields the stage efficiency, permitting the pressure ratio to be calculated. 
In a similar manner we can work backwards through the turbine stages. Thus 
we can calculate all the conditions on the far side of the L.P. rotor, which are 
then the entry conditions to the next stage, through which we can work in turn, 
and so on, integrating through all the compressor and turbine stages until the 
H.P. end is reached. 


Example 


To show the rdsults obtained by the method, Figs. 8 and 9 compare the 
measured and calculated performance, for an experimental turbocompressor 
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Fic. 9. Comparison of predicted and measured variation of compression temperature rise and 
efficiency with flame temperature. 


that was run at a constant value of .M, .p, and M = 1.0, while p was varied. 
p was varied by varying the turbine entry temperature, and the results are 
plotted against this variable. It is evident that the compression ratio is in 
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close agreement. The observed temperature rise was greater than that pre- 
dicted, which is the reason why the efficiency is low. The explanation of this 
can be seen by reference to the illustration of Fig. 1. Any leakage between 
stages from turbine to compressor will raise the apparent temperature rise. 
Considering this effect, the agreement is regarded as satisfactory. Fig. 10 
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Fic. 10. Comparison of predicted and measured value of compressor blade incidence. 


shows the calculated and observed values of the incidence on the compressor 
blades for all the rotors. In this figure two sets of predicted and observed 
results are shown, as predictions were made of the estimated results following 
a modification in the turbine blades to alter their angles. Again the agreement 
is fairly satisfactory. It is interesting to note the peculiarities of the observed 
results at rotors Nos. 10 to 14. It will be remembered that the theory used 
neglects any effect of stator angle. In other words it assumes that the stator 
angle can be always adjusted to suit operating conditions. In fact, of course, 
a compromise is necessary, and the lower curves show the effect of the error 
introduced. 


Application to Other Types of Unit 


The method described applies particularly to a certain type of turbocom- 
pressor unit. With suitable modifications it is applicable to any type of 
compound engine, having one or more compressor stages each of which is 
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independently driven by its own turbine stage. In such cases, it often happens, 
if there are relatively few stages, that the turbine stages may be operating 
under ‘‘choked”’ conditions which makes calculation easier, by controlling the 
mass flow ;M. 
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NOTES 


Removal of Traces of Oxygen from Gases 


For anaerobic studies in the Warburg apparatus, and for other purposes, it 
is often necessary to remove oxygen from commercial nitrogen or nitrogen— 
carbon dioxide gas mixtures. The method most commonly used consists in 
passing the gas over heated copper fragments in a combustion tube. The 
apparatus shown in the accompanying diagram and described below, has been 
found very convenient for this purpose. It is easily constructed, clean, rela- 
tively cool, easy to use, and occupies little laboratory space. Such apparatus 
has been used for several years in three laboratories. 


ASBESTOS PAPER 
Aspestos Pup 
ASBESTOS CoRD 


ASBESTOS PAPER 
WIRE 
ASBESTOS PAPER 





Fie. 1. 


A tube of Pyrex No. 172 glass about 70 cm. long, 16 mm. inside, 20 mm. 
outside diameter, is drawn down at one end to a tip suitable for connection to 
rubber tubing. The tube is well packed with copper oxide wire fragments, 
which are held in place by rolls of copper gauze. The greater part of the length 
of the tube is wrapped with two layers of thin asbestos paper and then wound 
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evenly, 3 to 3.5 turns per cm., with nearly 40 ft. of Nichrome wire, 2 ohms per 
ft. Two more layers of asbestos paper are wrapped on, and over this a layer 
of asbestos cord is wound closely. An asbestos paste is prepared by crumbling 
ordinary steam pipe lagging with water. This paste is smeared and moulded 
onto the tube over the asbestos cord to a thickness of about 15 mm. After 
this has dried it is protected by a further layer of asbestos paper tied on with 
wire. 

It is useful to make one or two windows in the windings so that the condition 
of the copper in the tube may be seen. To make these, one or two turns of the 
Nichrome wire and asbestos cord are spaced widely, and the layers of asbestos 
paper are cut away with a razor blade, except under the wire, before applying 
the asbestos paste to the remainder of the tube. The apparatus as described 
is used with 110 v. d-c. power source but 110 v. a.c. should be equally satis- 
factory. No external resistance is used except initially when it is advisable 
to put a lamp in series until the asbestos lagging has thoroughly dried. By 
means of a high temperature thermometer embedded in the copper fragments 
it was found that the temperature reached about 450°C. within 20 min. of 
turning on the current, and the maximum temperature was about 650° to 


700°C. 


To reduce the copper, hydrogen mixed with a carrier stream of nitrogen or 
other inert gas is passed through the hot tube. The tube takes an hour or 
more to cool after the power is disconnected and it is advisable to maintain a 
slight flow of inert gas through the tube until it is cool. 
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